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Statistical fluctuations in the transverse distribution of sources in rela- 
tivistic heavy ion collisions and an asymmetric emission profile associated 
with the wounded nucleons lead to rapidity dependence of the reaction 
plane. The size of this effect is estimated for the gold-gold collisions at 
the highest RHIC energy (y / ijv r /v = 200 GeV) in the Glauber model using 
GLISSANDO. The hydrodynamical evolution of the resulting torqued fireball 
is carried out in 3 + 1 dimensions considering a perfect fluid. Hadroniza- 
tion is simulated using THERMINATOR including non-flow contribution com- 
ing from resonance decays. Some experimental measures that can be used 
to detect the torque effect are proposed. 
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Fig. 1. (Color online) The inclusion of a random asymmetric cluster, schematically 
depicted in Fig. 1(a) twists the principal axis. Due to the asymmetry in the 



emission profiles of forward(backward) moving wounded nucleons, f±, seen in Fig. 



1(b) this twist is rapidity dependent (/ refers to the emission profile for binary 



collisions). The schematic figure of the torqucd fireball, Fig. 1(c) elongated along 
the m axis. The direction of the principal axes in the transverse plane rotates as 
7711 increases. The left and right pictures correspond to the rank-2 (elliptic) and 
rank-3 (triangular) cases, respectively. The effect occurs event-by-event. 



1. Introduction 

Dynamical properties of the early stage of relativistic heavy ion colli- 
sions can be probed using long range rapidity correlations. Starting from 
the wounded nucleon [H [2] approach to the initial stage of heavy ion colli- 
sions a torqued fireball picture naturally arises as a result [3]: the principal 
axes are twisted in opposite directions in the forward and backward regions. 
The underlying reasons for this are the statistical fluctuations in the trans- 
verse density of the sources (wounded nucleons) and the assymetry of the 
particle emission function (forward (backward) moving wounded nucl eons 



emit particles preferably in the forward (backward) direction). In Fig. 1(a) 
we represent schematically a cluster of wounded nucleons causing a twist m 
the principal axes. If we are dealing with a different number of forward and 
backward moving wounded nucleons, due to the asymmetry of the emission 
functions, this twist is rapidity dependent. The initial energy density is: 

F(v\\,x,y) = (1 - a)[p+(x,y)f+iv\\) + P-(x,y)f-(V\\)) + apbw{x,y)f(r)\\), (1) 

where p±, pbin are the density of forward (backward) moving wounded 
nucleons and binary sources. The respective emission profiles f±, f are 



depicted in Fig. 1(b) (in the chosen parametrization / = /+ + /- and is 
symmetric in rapidity). The parameter a corresponds to the fraction of 
binary sources in this mixed model [H [6] . 

The detection of this effect from the data consisting of the final momenta 
of observed particles in an event is a challenge as one still has to take into 
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account the effect of the fireball expansion and hadronization. 



2. Simulation results 

To evaluate the size of the torque effect in the initial fireball we perform 
a Monte-Carlo simulation using GLISSANDD [7J in the Glauber mixed model. 
In a given event, the angle of the principal axis for the Fourier moment of 
rank-/c for n sources at some 7711 is given by: 



*W = arctan 1 tttt ■ (2) 

Here the weight of the sources are: wi = (1 — a)/± for forward (backward) 
moving wounded nucleons and Wi = af for binary. We then define the dif- 
ferences: A FB = ^ (2) (f?||) - ^ (2) (-V\\) and A FB(BC) = ¥ 2 \±r] ll ) - ¥ 2 \0). 
The magnitude and sign of the resulting torque angle fluctuate from event 
to event. For Au+Au collisions at the highest available energy at RHIC 
the expected size o£ Ap B is of 20° for central and 10° for mid-peripheral 



collisions (see Fig. 2(b) ) 




Fig. 2. Distribution of the forward-backward angle difference for the elliptic de- 
formation, Aps, for the 20 — 30% centrality class in Fig. 2(a) for 7711 = 0.5 and 
7711 = 2.5 (the latter is the wider distribution). Rapidity dependency of the rms of 



Afb for several centrality classes Fig. 2(b) 



The 3 + 1 dimensional (perfect fluid) hydrodynamical simulation of the 
fireball expansion was done as in [8]. The results suggest that the effect 
surv ives this phase with almost no change in its size, as can be seen in Fig. 
3(a) The energy density used as a starting point for the hydrodynamical 
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simulation is given by: 



e(v\\,%,y) =(! - ^)[p+{Rx,Ry)f+{v\\) + P-{R T x,R T y)f-{-q\\)] + apbm(x,y)f(v\\), 

(3) 

where the operators R and R T rotate the density of the forward and back- 
ward moving wounded nucleons (respectively) by a fixed angle in opposite 
directions (5° for the centrality c = 20 — 25%). 

However, this effect is significantly washed away in the hadronization 
stage. Using the result from Appendix C of [9] we can estimate the deviation 
of the event angle from the fireball angle, as seen in Fig. |3(b)| Higher 
multiplicity as well as higher flow coefficient lower the statistical noise. 
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Fig. 3. Angle as a function of space-time rapidity for different values of proper time 
given by 3+1 perfect hydrodynamic evolution of the fireball can be seen in Fig. 3(a) 



Estimate of the difference between the fireball angle and the event angle for three 
different multiplicities and two different values of the flow coefficient (red, green 
and blue correspond respectively to n = 600, 100 and 20; solid lines correspond to 
Vk = 5% and dashed to 2%). 



3. Proposed measures: cumulants 



We propose experimental measures, which should be possible to apply 
with the available RHIC data, based on the use of cumulants involving 
particles from different rapidity regions. If confirmed, this would support 
the initial assumptions concerning the fluctuations in the early stage of the 
fireball formation and asymmetric rapidity shape of the emission functions 
of the moving sources in the nucleus-nucleus collisions. Two- and four- 
particle cumulants using particles from forward and backward rapidity bins 
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could be used for this purpose: 
cos(fcA FB ) {2} = X ' = = (cos(fcA FB )) events + nflow 



l e ik[(4>F,l+4>F,2)-(<t>B,l+4>B,2)]\ 

co S (2kA FB ) {4} = | eifc[( ^ 1 _^ a) _ ( ^ -1 _^ t2) / ) = (cos(2fcA i , s )) cvonts + nflow. 



(4) 



Another possibility would be to use particles from three different rapidity 
bins (forward, central and backward): 

/ e i2[(<f>F-<j>c,i)-(<t>B-<l>c,2)]\ _ / e i2[(<t>F-4>a,i)M<t>B-^>o,a)]\ 

AfbcH} = 9 

V2,FV2,BV^ C 

= (2sin(2Ai?c) sin(2A sc -)) cvcnts + nflow ps 8cov (A F cA B c) + nflow. 

(5) 

The results of the THERMINATOR simulation can be seen in Fig. and 
suggest that it should be possible to detect the effect of the torque. Full 
lines correspond to the expected result (without non-flow). 



4. Conclusions 

We have proposed experimental measures that should be able to detect 
a torqued fireball effect resulting from the fluctuation of the sources trans- 
verse position and their asymmetric emission profile. Due to the decrease in 
statistical noise with higher statistics and flow coefficient this effect should 
be looked for in mid-central or mid-peripheral collisions. Since the effect in- 
creases with rj, the experimental search at RHIC should involve the forward 
and backward TPCs. 
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Fig. 4. Pseudo-rapidity dependence of 2-particle (top row) and 4-particle cumulant 
using forward-backward rapidity bins (middle row) and using three rapidities (bot- 
tom row). In the left panels only primordial particles are considered whereas in the 
ones on the right resonance decays are included. In the latter case only charged 
pions, kaons, protons and antiprotons are considered and a cutoff in transverse 
momentum is included (450 MeV < px < 3 GeV). 



[3] P. Bozek, W. Broniowski and J. Moreira, Phys. Rev. C83 (2011) 034911, 
1011.3354. 

[4] D. Kharzeev and M. Nardi, Phys. Lett. B507 (2001) 121, nucl-th/0012025. 
[5] PHOBOS, B.B. Back et al., Phys. Rev. C65 (2002) 031901, nucl-ex/0105011. 
[6] PHOBOS, B.B. Back et al., Phys. Rev. C70 (2004) 021902, nucl-ex/0405027. 
[7] W. Broniowski, M. Rybczyhski and P. Bozek, Comput. Phys. Commun. 180 
(2009) 69, 0710.5731. 

[8] P. Bozek and I. Wyskiel, Phys. Rev. C79 (2009) 044916, 0902.4121. 



sqm2011'jmoreira printed on March 30, 2012 



[9] W. Broniowski, P. Bozek and M. Rybczyriski, Phys. Rev. C76 (2007) 054905, 
0706.4266. 



